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SPHERICAL  DIFFUSION  FLAMES 


Lav  has  shown  that  the  analysis  of  spherical  diffusion  flanes  is  quite 
similar  to  that  of  counterflov  diffusion  flames,  so  that  some  explanation 
is  needed  for  devoting  a  separate  lecture  to  them.  There  axe  tvo  good 
reasons.  First,  the  constant-density  approximation  has  been  used  through¬ 
out  these  lectures  in  discussing  all  but  plane  flames,  so  there  is  room 
for  a  problem  which  does  not  neglect  variations  in  density.  (Plane 
diffusion  flames  have  to  be  chambered,  i.e.  the  reactants  must  be  supplied 
at  finite  locations,  which  leads  to  distracting  complications.)  Secondly, 
the  spherical  diffusion  flame  can  lead  to  quite  different  (and  unusual) 
responses.  These  arise  in  the  technologically  important  application 
to  the  quasi-steady  phase  of  fuel-drop  burning,  when  a  more  realistic 
boundary  condition  than  the  conventional  one  is  used. 

1.  3asic  Equations 

The  fuel  is  supposed  to  be  supplied  as  a  liquid  at  the  surface  of  a 

sphere  of  radius  a  (figure  1),  where  the  heat  from  a  concentric  flame 

sheet  evaporates  it.  In  turn,  the  flame  sheet  is  sustained  by  the 

reaction  of  the  gaseous  fuel  and  oxidant  in  the  ambient  atmosphere.  If 

the  representative  mass  flux  M  is  taken  to  be  X/c  a,  then  a  is  the 

r  P 

length  unit.  The  density  at  infinity  will  be  used  for  p^. 

We  seek  a  spherically  symmetric,  steady  solution  of  the  full  equations 
(2.18-20)  under  the  assumption  of  unit  Lewis  numbers  for  the  reactants. 

The  equation  of  continuity  integrates  immediately  to 

2 

r  pv  =  M  (const.),  (l) 

where  !!  is  loosely  called  the  burning  rate.  (The  liquid  fuel  evaporates 


-9.2- 


at  the  rate  but  the  gaseous  fuel  only  burns  at  that  rate  if  there  is 

no  outflow  of  fuel  at  infinity.)  Although  the  radial  velocity  v  is  not 
determined  by  this  result,  the  thermal-chemical  equations  (2.20)  are 
uncoupled  from  the  fluid-mechanical  equations  (2.l8b,19)  because,  under 
steady  conditions,  they  only  involve  the  mass  fluxpv,  which  is  determined 
(as  M/r  ).  Once  T  is  found,  the  density  follows  from 


pT  =  T 

CO 


(2) 


and  then,  v  is  known;  the  momentum  equation  (2.19)  serves  only  to  determine 
the  pressure  field. 

We  have  to  deal  once  more  with  equations  (8.2),  where  now,  however. 


_M_  _d_  _1 _ d_  _d_ 

2  dr  2  dr  '  dr  ’ 
r  r 


(3) 


The  boundary  conditions  are 

T  =  T  ,  dT/dr  =  ML,  X-M_1dX/dr  =  0,  Y-M^dY/dr  =  1  at  r  =  1,  (k) 


T  +  T  ,  X  X  ,  Y  -►  0  as  r  «,  (5) 

00  00 

where  Tg,  L,  T^,  X^  are  given.  Prescription  of  the  surface  temperature 

7_  will  later  be  replaced  by  the  requirement  of  liquid-vapor  equilibrium; 

the  latent  heat  of  evaporation  L  is  positive;  the  conditions  (kc,d)  ensure 

that  the  sphere  is  a  source  of  fuel,  but  neither  a  source  nor  sink  of 

oxidant  (or  anything  else);  the  prescribed  oxidant  fraction  at  infinity 

rust  of  course  satisfy  0  <  X  <1;  and  the  condition  (5c)  ensures  that 

00  ■” 

all  the  fuel  has  originated  at  the  supply. 

The  sixth-order  system  of  differential  equations  is  therefore  subject 
to  seven  boundary  conditions,  as  a  consequence  of  which  M  is  determined 
as  a  function  of  V  (depending  also  on  T_,  L,  T  ,  and  X  ).  The  maximum 
temperature  could  be  used  to  characterize  the  solution,  as  for  the  counterflow 


diffusion  flame,  but  ve  shall  use  M  instead.  The  burning  rate  is  of 
greater  interest  in  applications,  and  it  arises  naturally  in  the  analysis. 
The  Shvab-Zeldovich  relations  are 

T+2X  =  (T  ■ -L)(l-e"M/r)+(T  +2X  )e_l^r  =  Z  (M),  (6) 

5  oo  op  • 

T+2Y  =  T  +(T  -T  )(l-e“M^r)  =  Z(M),  (7) 

CD  a.  CO  +  t 

where 

T  =  T  -L+2 
a  s 

is  called  the  adiabatic  flame  temperature.  For  =  T  ,  no  heat  is 
conducted  to  or  from  the  ambient  atmosphere  when  (as  is  usually  the  case) 

Y  =  o(r-1)  as  r  ->  «;  this  conclusion  follows  from  the  relation  (T), 
which  gives 

lim(Uirr2dT/dr)  »  lim  [lnrr2(T  -T  )  ( -Me_M/r/r2 )  ]  =  UirM(T  -T  ) .  (8) 

a  oo  w  a. 

r-*» 

Me  further  see  that 

T  >  T  (9) 

a 

will  ensure  the  conduction  of  heat  to  the  environment,  the  aim  of  combustion 
in  practice. 

The  problem  now  reduces  to  one  for  T  alone,  namely 

1(7)  =  PXYe-9/*,  T  =  T  ,  dT/dr  =  ML  at  r  =  1,  T  +  as  r  +  ^,(10) 

where  X  and  Y  are  to  be  suppressed  in  favor  of  T  by  means  of  the 
Shvab-Zeldovich  relations  (6),  (7).  We  shall  first  derive  the  frozen 
■*  0)  and  equilibrium  (Pe“9|/T  ->  ®)  limits  of  the  solution. 


?or  frozen  combustion,  start  by  setting 


e  =  IPX  e 


ecual  to  zero.  We  find 


M  (l-l/r) 

=  T  -L+Le  w  vith  M  =  in[l+(T  -T  )/L] 

5  W  oo  s 


..’here ,  since  Y  must  be  positive. 


M  >  0,  i.e.  T  >  T  . 
w  ’  oo  s 


This  is  the  pure  evaporation  solution  vith  no  combustion;  the  environment 
must  be  hotter  than  the  surface  in  order  to  supply  the  heat  necessary  to 
.-aporize  the  liquid.  (There  is  no  condensation  solution.)  The  result  is 
:ot  uniformly  valid,  since  for  r  =  0(e_1)  the  reaction  terra  is  comparable 
:o  the  convection  terra  in  the  equation  (10a).  The  variable 


R  =  er 


.eads  to  the  expansion 


t  -  v^w[vvL-2(1-e  )1/R+*-- 


7r.s  frozen  limit  is  not  an  extinguished  state,  but  rather  one  in  which  all 
he  reaction  takes  place  at  essentially  constant  temperature  far  from  the 
upply . 

Turning  r.ov  to  the  (Burke -Schumann)  equilibrium  limit ,  we  use  the  Shvab-ZeldovicJ 
elations  (6),  (7)  to  obtain 


'ZJ'AJ 

°  1 

T  ». 

1 

V 

II 

-M  / r  »  Y  =  ( 

Z  (>'  ) 

,  +  e 

[(!«.)•  6  -1  ( 

-M  /r 

1-(1+X  )e  e 


for  r  <  r*,  (l6) 


here  the  boundary  condition  (10b)  fixes 


7*  =  T  +  yX  with  y  =  (T  -T )/(l+X ).  (19) 

"  00  00  a  00  OO 

For  consistency* the  value  (l8)  of  r*  must  be  greater  than  1,  i.e. 

T  -7  >{L-2)/X  .  (20) 

oo  g  oo 

Thiie  this  condition  is  automatically  met  when  the  inequalities  (9),  (13b) 
are  satisfied  (as  whey  are  ir.  the  practically  important  case  to  which  we 
shall  limit  our  later  discussion)  it  is  of  interest  to  note  that  there 
is  a(  second  limit,  found  by  Buckmaster,  when  the  condition  is  violated. 

The  flame  sheet  lies  at  the  surface  itself,  instead  of  some  distance  away, 
so  chat  there  is  no  oxidant-free  region. 

The  flame-sheen  structure  in  the  Burke-Schumann  limit  is  the  same 
irrespective  of  whether  the  count erf lev  or  spherical  diffusion  flame  is 
being  considered  (cf.  end  of  section  8.1).  The  continuity  conditions 

5(7)  =  c  ( M )  =  6(Y)  =  0  (21) 

are  required  for  a  structure  to  exist,  but  the  gradients  of  T,  X>  and  Y 
aru  different  on  the  two  sides  of  the  flame  sheet  because  the  latter  is  a 
source  of  heat  and  a  sink  of  both  reactants.  The  Shvab-Zeldovich  relations 
show,  however,  that 


5{dT/dr+2dX/ir)  *  <5  (dT/dr+2dY/dr )  =  0. 


(22) 


-►  oo 


We  turn  no w  to  the  full  response  curve  M( V)  in  the  limit  9 
: ir.g  once  again  that  the  burning  rate  is  a  more  convenient  and 
naif  icanc  parameter  than  the  maximum  temperature  to  characterize  the 
.ution.  Here  the  response  is  S-shaoed  if  T  ,  L,  T  satisfy  the 
•qualities  (9)  and  (13b),  i.e.  when  the  heat  flux  (8)  is  zo  the 
•ironment ,  as  required  in  practice,  and  the  surface  temperature  is  less 
m  that  of  the  ambient  atmosphere.  The  transition  from  an  S-shaped 
jpor.se  to  a  monotonic  one  occurs,  for  large  activation  energy,  when 
;  temperature  gradient  (i.e.  the  heat  flux)  beyond  the  flsuae  is  small, 
that  it  nay  be  described  by  setting 

T  -T  =  k/e  with  k=  const.;  (23) 

&  00 

■e  T  and  L  are  supposed  fixed  as  T  varies.  The  inequality  (13b) 
,ch  ensures  a  weak-buming  branch  for  the  response,  requires 

i  <  2.  (2h) 

e  that  the  requirement  (23)  for  the  Burke-Schumann  limit  is 
omatically  satisfied,  ensuring  a  strong-burning  branch  of  that  form. 

Equilibrium  will  be  assumed  behind  the  flame  sheet  even  though  the 
per  at 'ore  does  not  rise  above  Ttt  (=  T^,  to  leading  order)  there  by 
3(l)  amount.  The  combustion  is  frozen  between  the  surface  and  the 
me  sheet,  because  of  the  requii'ment  (13b).  The  Shvab-Zeldovich 
ation  (7)  therefore  gives  us 

T  =  T3j+(k/9)(l-e'M/r)  for  r  >  r*,  (25) 

-/lew  of  the  assumption  (23);  and  the  boundary  conditons  (lOb,c)  lead  to 


for  r  <  r. 


(26) 


L 


-9.7- 

WI(1-1/r) 


#* 


These  formulas  are  correct  to  any  order  in  9  ^ ,  provided  that  M  is 
determined  to  the  sar.e  order;  we  shall  only  need  a  leading-order  determine 
tion  of  M,  and  that  will  be  understood  in  what  follows.  Leading-order 
continuity  of  T  now  shows  that 


r*  =  M/[M+tn(L/2}],  (2?) 

a  result  that  exhibits  the  need  for  the  inequality  (2k).  To  analy?  he 
reaction  zone  we  shall  also  need  the  leading-order  result 

X*  =  (l+Xje-M/r*-l,  (23) 

coning  from  the  Shvab-Zeldovich  relation  (6). 

Determination  of  M  comes  from  analysis  of  the  reaction  zone,  for 
which  the  appropriate  variable  is 

5  =  e(r-r*).  (29) 


Coefficients  in  the  layer  expansion 

T  =  T  -0_1T2(*>  +...  with  6  =  (l/T)..  (30) 

ao  oo  _L 

are  considered  to  he  functions  of  £•  The  Shvab-Zeldovich  relation  gives 

,  ,  -M/r»  2 

Y  =  T~(p-«*)+ —  with  =  -k(l-e  J/T^,  (31) 

so  that  the  structure  is  governed  by  equation  (8.2U)  with  Tf+  replaced 
by  T  .  As  there,  the  structure  equation  provides  an  expression  for  the 
■enperature  gradient  in  the  field  (26)  in  front  of  the  flame  sheet,  whence 


find  the  equation 


T  /KC,  e 
00  “ 


for  M  as  a  function  of  P.  This  function  can  be  defined  parametrically, 
as  for  the  counterflow  diffusion  flame  at  the  end  of  section  8.2.  Thus, 

M  is  given  as  a  function  of  r*  by  the  formula  (27),  and  then  V  is 
determined  as  a  function  of  r#  by  equation  (32)  with  the  substitutions 
(28)  and  (31b)  for  X*  and  $#• 

The  corresponding  response  curve  is  shown  in  figure  2  for  several 
values  of  k.  The  curve  joins  the  frozen  limit  (12b)  to  the  Burke -Schumann 
limit  (17),  monotonically  if  k  is  not  too  positive  but  otherwise  via  an 
3.  Reponses  in  the  shape  of  an  S  appear  to  be  associated  with  a  flux  of 
heat  from  the  flame  sheet  to  an  environment  that  is  hotter  than  the  supply 
surface,  the  normal  state  of  affairs  in  practice.  However,  it  has  never 


been  proved  that  the  inequalities  (9)  and  (13b)  will  ensure  an  S-response 
in  the  limit  0  •>  «  . 


General  Extinction  and  Iznition  Analyses 


T 


also , involving  the  perturbation  M^.  It  provides  the  stronger  matching 
condition 


b  =  -£-M,+.. .  as  with  =  (l+y)?-^/?* 


(35) 


in  the  canonical  problem  (3.32,33),  at  which  the  analysis  finally  arrives 


Cl 


-  j,  a  p  -s 

?  =  1-Y,  Ve  =  r*  T°  PAtlVe  ; 


(36) 


here  r#,  T^,  and  Y  have  the  definitions  (l8)  and  (19). 

Since  the  problem  is  well  posed  under  the  weaker  condition,  it  follows 

that 

M.  =  -lin  (J+i)  (37) 

car.  be  calculated  as  a  function  of  V  ,  thereby  determining  the  response. 

e 

This  response  is  shown  in  figure  3  for  various  values  of  Y»  it  has  been 
found  numerically,  but  not  yet  been  proved,  that  the  curve  turns  to  form 
a  C  whenever  y  is  positive,  consistent  with  the  conjecture  that  the 
inequality  (?)  must  hold  for  an  S-response. 

The  general  ignition  analysis  that  was  promised  in  section  8.6  will 
r.ow  be  presented.  The  starting  point  is  the  analog  of  the  assumption 
(33),  namely 

M  =  M  +  e“1i41+. .. ,  (38) 

where  M  has  the  frozen  value  (12b).  Correspondingly,  the  reaction  all 

V 

takes  dace  far  from  the  surface,  the  combustion  field  being  frozen  to 


all  orders  at  any  finite  r.  (The  latter  is  ensured  by  the  inequality  (13b).) 


.he  formulas 


"-  =  T  -L+Le 
0  s 


M  (1-1/- ) 


,  Tx  =  ^Ld-l/rle 


M  (1-1/r) 


(39) 


re  obtained  by  satisfying  conditions  (lOb,c)  at  the  surface;  but  these  do 
ot  satisfy  the  condition  (iOd)  at  infinity. 

The  appropriate  variable  at  large  distances  is 


R  =  r/9. 


(1*0) 


id  then  coefficients  in  the  expansion  (30)  are  to  be  considered  functions 
'  R.  The  equation  for  <j>  is 

-  ]e  0  with  P  =  iPX  9  e  ; 

.  V  OO 

T  R 


(T  -T  )M 

h  s  <5‘  ■>+ 


(1.1 ) 


is  to  be  solved  under  the  boundary  conditions 

y* 

0  =  Le  V'(-M-  +M  /R)T2  +...  as  R  -►  0,  $  *=  o(l)  as  H>», 

i.  tt  00 


(1*2) 


.e  first  of  which  cones  from  matching  with  the  expansion  coefficients 
?)  for  finite  r. 

Since  the  problem  is  well  posed  under  the  weaker  condition 

M 


lin  R$  =  Le  VM  /T2, 

„  CO 

R-*0 


(1*3) 


follows  that 


?  M 

.v  =  lim  (M  /R  -  T“d/Le  w) 
T  R-0  w 


(1*4) 


n  be  calculated  as  a  function  of  P  ,  thereby  determining  the  response. 

^  Mw  2  2  2  ^1 

e  relation  between  the  scaled  parameters  (Le  /T  )M  ,  (M  L  e  w/t  )P 

00  J.  V  oo  V 
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(^5) 


figure  4  shows  the  corresponding  response  for  several  positive  values  of  8. 
It  is  found  numerically ,  but  has  never  been  proved,  that  the  curve  does  not 
uurr.  for  non-positive  values  of  8,  consistent  vit'n  the  conjecture  that 


the  ir.ecuality  (9)  must  hold  for  an  3-resoonse. 

3  e/(Tanxo) 

It  is  worth  noting  that  V  is  0(8  e  °°  )  at  extinction 

6/T 


.-2 


compared  to  0(0  e  )  at  ignition,  so  that  I  lies  to  the  left  of  I, 

as  required.  The  part  between  the  ignition  and  extinction  points,  i.e.  the 
middle  branch  of  the  S,  will  not  be  discussed  here  because  it  is  very  similar 
to  that  for  the  counterflow  diffusion  flame  (section  8A.).  Again  the 
partial -burning  and  premixed  flames  occur,  with  the  same  structures  for 
their  reaction  zones.  Since  the  stability  analysis  of  the  branch  is 
confined  to  the  reaction  zones,  it  is  the  same  here  as  in  section  8.5, 
so  no  further  discussion  is  necessary. 


4 .  Surface  equilibrium 

There  are  at  least  two  ways  in  which  the  fuel  can  be  supplied  as  a 
liquid  at  the  surface:  the  sphere  can  be  completely  liquid,  i.e.  a  fuel 
drop,  or  it  can  be  a  liquid-saturated  porous  solid.  Whatever  the  method  of 
supply,  it  is  difficult  to  justify  prescription  of  the  surface  temperature 
!_  unless  we  abandon  the  specification  of  L.  Maintaining  a  value  of  T 
would,  in  general,  require  heating  or  cooling  the  liquid  at  the  surface, 
and  that  would  upset  the  heat  balance  represented  by  the  boundary  condition 
(13c).  Prescription  of  surface  temperature  should,  more  realistically,  be 
reclaced  by  the  requirement  of  liquid-vapor  equilibrium  at  the  surface. 


i.e.  the  Clausius-Clapeyron  relation 


Y 


(T./Tj 


B  S/V5/Ts 

e 


with  8,0  consts..  (U6) 


b 


Here  is  the  boiling  temperature,  i.e.  the  value  of  Tg  for  which 


Y  =  1;  it  is  related  to  the  pressure  level  in  the  combustion  field  by 

S  "§/Tb 

T?e  =  p  /k  with  k  =  const., 

b  c 


(U7) 


In  order  to  transform  the  Clausius-Clapeyron  relation  into  a  temperature 

condition,  ve  eliminate  Y  between  it  and  the  Shvab-Zeldovich  relation 

s  » 

(7)  to  obtain 


.  e/T, -e/T 

2<VV  e 


=  2-L+(T  -T  +L-2)e 

CO  s 


-M 


(U8) 


This  replaces  the  boundary  condition  (10b' ;  clearly  Tg  changes  with  M. 

Determination  of  the  response  M(V )  now  requires  that  T  be 

s 

calculated  afresh  for  each  point.  In  general,  iterations  are  involved 

because  the  structure  problem  determing  M  as  a  function  of  V  contains 

Is  (see,  for  example,  the  definitions  (36)).  Different  results  are 

obtained  according  to  the  way  in  which  V  is  varied,  via  pc  or  the  radius 

a  of  the  sphere,  the  reason  being  that  pc  appears  in  the  boundary 

condition  (1*8)  through  T^.  (Since  pr  has  been  taken  proportional  to 

n  ,  and  M  inversely  nronortional  to  the  radius  a  of  the  sphere,  the 
-  c  r 

2 

definition  (2.23b)  shows  that  V  is  proportional  to  p^a  .) 

Normandia  and  Ludford  have  considered  the  response  when  V  is  varied 
via  a.  In  particular,  they  find  that  the  curve  is  again  S-shaped  when  the 
inequalities  (9),  (13b)  are  satisfied  by  the  Tg  corresponding  to  small  V, 
i.e.  the  solution  of  equation  (1*8)  when  M  has  the  value  (12b).  (There 
is  always  such  a  solution,  satisfying  the  left  inequality  automatically.) 

The  previous  analysis  (for  Tg  fixed)  is  qualitatively,  but  not 
quantitatively,  correct:  while  extinction  and  ignition  analyses  can  be 
performed  as  in  section  3,  figures  3  and  •*  cannot  be  carried  over  because 
the  variations  in  T  .though  only  0(0  1)  on  the  upper  and  lower  branches 


of  the  S-res?onse,  modify  the  extinction  and  ignition  values.  These 
modifications  have  been  worked  out  by  Normanaia  and  Ludford. 

The  results  are  strikingly  different  when  V  is  varied  via  ?c> 
as  Janssen  &  Ludford  (1983a)  have  shown.  The  S  is  replaced  by  the 
rather  odd  shapes  in  figure  5,  which  covers  the  practical  cases  of 
heat  being  conducted  to  the  environment.  These  responses  -were  so 
unexpected  that,  to  obtain  more  confidence  in  their  validity,  numerical 
integrations  of  the  problem  (I0a,c,d,  18)  were  performed  for 
3  =  10.  In  every  instance,  the  numerical  results  confirmed  the  essential 
features  of  the  asymptotic  responses .  Nothing  similar  has  apparently 
been  obtained  in  combustion  theory,  and  certainly  not  in  previous 
studies  of  diffusion  flames. 

Two  features  of  these  responses  deserve  to  be  pointed  out ,  since 
they  contradict  conventional  wisdom.  First,  the  3urke-Schumann  value 
(IT)  is  not  attained  on  any  of  the  curves  as  V  »  .  The  3urke- 
Schumann  solution  discussed  in  section  1  is  considered  a  good  approxima¬ 
tion  since  in  practice  0  is  large.  3ut,  instead  of  standing  off 
from  the  sphere  for  high  pressures,  separating  equilibrium  regions, 
the  flame  sheet  actually  moves  to  the  surface,  forming  the  second 
(Buckmaster)  equilibrium  limit  mentioned  in  connection  with  the 
condition  (20).  Secondly,  the  burning  rate  decreases  over  most  of 
each  response,  whereas  the  general  belief  is  that  it  should  increase. 
Negative  slope  of  the  response  curve  is  thought  to  be  inevitably 
associated  with  instability,  and  indeed  we  found  that  to  be  the  case 
for  the  S-responss  at  the  end  of  section  3.  But  experience  refutes 
such  a  conclusion  here:  by  and  large,  fuel  drops  do  burn  steadily. 


A  physically  reasonable  range  of  pressures  focuses  attention  on 
the  left  ends  of  the  curves  in  figure  5>  and  there  the  decrease  of 
burning  rate  uith  pressure  is  by  no  means  insignificant.  Moreover, 
the  computed  response  curve  for  net hand  over  the  range  .002  to  500 
at- o spheres  clearly  shovs  the  steady  decline  (Janssen  &  Ludford  1983b). 
Experiment  should  therefore  be  able  to  give  a  clear-cut  decision  on 
the  physical  reality  of  this  'unexpected  phenomenon. 
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Figure  Captions 


9.1  Burning  fuel  drop. 

9.2  Steady-state  responses  when  the  combustion  is  nearly 
adiabatic.  Drawn  for  L  =  0.2  ,  T  -  0.2  ,  X  «  1  . 

00  09 

9.3  Extinction  curves..  For  y  ^  0  no  turning  point  is  found 
numerically. 

9.4  Ignition  curves.  For  B  _<  0  no  turning  point  is  found 
numerically. 


9.5  Sketch  of  steady  state  responses  when  there  is  surface 
equilibrium. 
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